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Charge stripes and spin correlations in copper-oxide superconductors 

J. M. Tranquada^* 

'^Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 

Recent neutron diffraction studies have yielded evidence that, in a particular cuprate family, holes doped into 
the Cu02 planes segregate into stripes that separate antiferromagnetic domains. Here it is shown that such a 
picture provides a quantitatively consistent interpretation of the spin fluctuations measured by neutron diffraction 
in Lai.85Sro.i5Cu04 and YBa2Cu306+a^. 



1. INTRODUCTION 

To obtain superconductivity in a layered 
copper-oxide compound, it is necessary to intro- 
duce charge carriers into the antiferromagnetic 
Cu02 planes. Recent neutron diffraction studies 
of the system Lai.6-2;Ndo.4Sra;Cu04 [|l],|j provide 
evidence that the dopant-induced holes choose to 
segregate into periodically-spaced stripes which 
separate antiferromagnetic domains, in a man- 
ner similar to that found in hole-doped La2Ni04 
|3|. The charge and spin stripe modulations are 
identified by the appearance of scattering at in- 
commensurate positions. In the Nd-doped sys- 
tem, elastic scattering is observed, correspond- 
ing to static stripes. In pure La2_a;Sra;Cu04, the 
magnetic scattering that is observed is purely in- 
elastic Q. Where samples with and without Nd, 
but with the same Sr concentration, have been 
measured, the incommensurate (IC) splittings of 
the magnetic signal are found to be essentially 
identical [HJ|J§1. It has been proposed that the 
spin correlations in the two systems are funda- 
mentally the same, thus implying similar charge 
correlations. The static nature of the stripes in 
the Nd-doped system is attributed m to pinning 
of the otherwise dynamic correlations by a spe- 
cial distortion of the lattice . That distortion 
is driven by purely ionic interactions and is stabi- 
lized by the smaller ionic radius of the substituted 
Nd. 

To strengthen the case for the charge-stripe in- 



terpretation, it is necessary go beyond a simple 
comparison of incommensurabilities. In this pa- 
per I will show that the neutron scattering mea- 
surements of spin fluctuations in superconducting 
Lai.85Sro.i5Cu04 and YBa2Cu306.6 are quanti- 
tatively consistent with what one should expect 
based on the stripe picture. The general expecta- 
tions are described in section 2. Specific analyses 
for the two different superconductors are given in 
sections 3 and 4, respectively. Section 5 contains 
some further discussion. 

2. MODEL 

First recall the situation for the 2D Heisenberg 
model 0. At T = 0, where long-range order ex- 
ists, the spin-wave dispersion is linear at small q 
(= Q — Qaf), and can be expressed as a; = cqQ. 
At finite temperatures, thermal fluctuations re- 
sult in a finite spin-spin correlation length ^. Spin 
waves are overdamped for q < 1/^, but spin- wave 
dispersion should still be detectable for g > 1/^ 
(as long as kT is less than the superexchange en- 
ergy J). The characteristic energy Fq separating 
the over- and underdamped regimes is then given 
by 



Fo w hco/^. 



(1) 
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As for the dynamical susceptibility, x"(Q,w), Fq 
acts somewhat like a gap energy Q] . If x" is inte- 
grated over Q to give x"{uj), then there is a peak 
at huj ~ Tq, with x"i^) going to zero as to goes 
to zero. 

Next consider an ordered stripe phase. The 
charge stripes need not affect the superexchange 



energy within the antiferromagnetic domains, but 
they certainly wiU weaken the magnetic coupling 
between domains. Castro Neto and Hone |g] 
have considered the effect of the anisotropic cou- 
pling on the spin wave velocity. For a large 
anisotropy, the velocity c for propagation paral- 
lel to the stripes is equal to co/V2; it is more 
difficult to estimate an effective velocity perpen- 
dicular to the stripes because the stripe modu- 
lation should result in several optical branches 
as well as the acoustic mode. An experimental 
test of this situation is given by recent measure- 
ments of spin waves in stripe-ordered La2Ni04.i33 
H. That study found c = O.Gcq for propagation 
parallel to the stripes. Dispersion was less well 
defined but comparable in the perpendicular di- 
rection. (Note that the spin waves disperse out of 
the IC points in reciprocal space that characterize 
the magnetic order.) Thus, it seems reasonable 
to take c — co/-\/2 as an estimate of the effective 
spin-wave velocity in an ordered stripe phase. 

Finally, consider a stripe phase with a finite 
spin-spin correlation length induced by fluctua- 
tions of the charge stripes. The characteristic en- 
ergy separating over- and underdamped modes 
should now be 



r « hco/V2t 



(2) 



Excitations at Hlo > F should look spin-wave-like. 
On the other hand, a scan at a fixed excitation 
energy smaller than F through the IC points char- 
acterizing the stripe correlations should show a 
broad peak centered at each IC point, with the 
width varying inversely as ^. To obtain an ana- 
lytic formula describing the g-dependence of the 
scattering, consider a model in which the spins 
have a uniform correlation within each domain, 
but where the correlations between domains fall 
off exponentially with distance. If the charge 
stripes are separated by n lattice spacings, then 
the intensity along Q = (/i, i,0) should be pro- 
portional to 



I{h) ^ \Ff 



l-p2 



(3) 



1 + p2 _ 2p cos(27rn/i) ' 

where p = (— 1)"+^ exp(— na/^), and F is the 
structure factor for a single domain. (The factor 
of —1 included in p takes care of the antiphase 



relationship between neighboring domains.) For 
the calculations below, it will be assumed that the 
spin correlations also fall off exponentially paral- 
lel to the stripes, with the same correlation length 

3. Lai 85Sro.i5Cu04 

Recently, Hayden et al. ||l^ reported inelastic 
neutron scattering measurements of x"(Q,a;) in 
Lai.86Sro.i4Cu04, covering a very large energy 
range. They found that x" a-t ?iw > 100 meV 
could be modelled reasonably well using spin 
wave formulas, with only a slight reduction in J. 
At lower energies, they showed that x"{^) has a 
peak at F = (22 ± 5) meV. If this latter result is 
combined with the value hcQ = (850 ± 30) meV- 
A |jll|, then by inverting Eq. (2) one obtains the 
estimate ^ = (27 ± 6) A = (7.0 ± 1.5)a. 

Figure H shows an inelastic scan, at a fixed en- 
ergy transfer of 3 meV, through the IC peaks 
found in a single crystal of Lai.85Sro.i5Cu04 Q. 
The peak positions are consistent with a charge 
stripe period of n = 4. The solid line through the 
data points is a calculation using Eq. (3) and the 
estimate ^ = 7a — only the amplitude and the lin- 
ear background have been adjusted to fit the data. 
(The calculated intensity, taking into account sig- 
nal from all 4 IC peaks, was also convolved with 
the experimental resolution function.) The calcu- 
lation clearly gives an excellent description of the 
data. 

4. YBa2Cu306.6 

In contrast to the 214 system, measurements of 
x" in YBa2Cu306.6 do not show well resolved IC 
peaks ||l^. As shown in Fig. 2, the scattered in- 
tensity is peaked at Qaf, but with a peak shape 
that has a rather flat top and steep sides. De- 
spite the lack of clear evidence for incommen- 
surability, it is possible nevertheless to describe 
the data with the same stripe model discussed 
above. To give the stripe picture an honest test, 
^ will again be estimated using Eq. (2). The 
value of F is approximately (30±5) meV Q, and 
hco = (670 ± 30) meV-A @,0, which together 
yield ^ = (16 ± 3) A = (4.1 ± 0.8)a. To estimate 



o 20 




YBa,Cu30,_^ T^ 



53 K 



150 



100 



50 




I I I I I I I I I I I I I I I 
2 0.3 0.4 0.5 
h 



0.6 



0.7 



Figure 1. Inelastic scan zXTilo — ?> meV through 
the incommensurate points (see inset) measured 
at T = 40K onacrystalof Lai.85Sro.i5Cu04 with 
Tc = 37 K. Sohd hne is the calculation described 
in the text; dashed line is the fitted background. 
Scan is along Q = (/i, i, 0), within the (hMS) zone, 
as indicated schematically in the inset. 



Figure 2. Inelastic scan at Tiw = 15 meV along the 
zone diagonal (see inset) measured at T = 10 K 
on a crystal of YBa2Cu306.6 with T^ = 53 K. The 
solid line is the calculation described in the text. 
The scan was measured along Q = (ft,, /i,— 1.8) 
within the (hhX) zone. The inset shows the scan 
direction within the (ft,, /c,— 1.8) zone. 



the stripe spacing, assuming that the hole den- 
sity within the charge stripes is the same as in 
the 214 system [n], it is first necessary evaluate 
the hole concentration within a plane. Inverting 
the formula ^M 

82.6(p-0.16)^ (4) 



-^ cj -^ c^r. 



= 1 



yields p « 0.088 for the Tc = 53-K sample, which 
corresponds closely to n = 6. The solid line 
through the data in Fig. 2 represents a calcula- 
tion using the same stripe model as above, with 
the parameter values n — Q and ^ = 4.1a. The 
calculated function was convolved with the spec- 
trometer resolution function, and only the am- 
plitude and background were adjusted to fit the 
data. 

For the calculated curve, the filling in of weight 
between the IC peaks occurs because ^ is less than 
the stripe spacing na. In this case the weak cor- 
relation between neighboring antiphase domains 



results in little interference. At Qaf the scat- 
tering from an individual domain does not get 
cancelled by the weak contributions from neigh- 
boring domains. In contrast, the situation for 
Lai.85Sro.i5Cu04 is ^ > na. That is the nec- 
essary for condition for obtaining well resolved 
IC peaks. Although the correlations between an- 
tiphase domains are perhaps less well defined in 
123, the smaller value of ^ suggests greater stripe 
fluctuations, a condition that appears to correlate 
with an increased Tc ||] . 

5. DISCUSSION 

It has been shown that the Q dependence of 
scattered intensity from low-energy spin fluctua- 
tions in Lai.85Sro.i5Cu04 and YBa2Cu306.6 can 
be consistently interpreted on the basis of dy- 
namic antiphase antiferromagnetic domains. The 
correlation lengths ^ used to model the scattering 



measurements were calculated from experimen- 
tal values of F and Hcq. It was pointed out that 
the IC peaks become unresolvable when ^ < na, 
where na is the charge-stripe spacing. 

There are, of course, other neutron scatter- 
ing results that are consistent with the stripe 
picture. A comparison of x"{uj) in super- 
conducting YBa2Cu306.6 and antiferromagnetic 
YBa2Cu306.i5 shows that, despite a difference in 
uj dependence, the overall weight within the range 
5 to 50 meV is comparable [H. Similar results 
were found in a comparison of Lai.85Sro.i5Cu04 
and La2Cu04 |10|, although the spectral weight 



for the superconductor is somewhat reduced when 
the comparison is extended up to 300 meV. With 
a strong segregation of holes of density nt, the 
density of Cu moments contributing to the low- 
energy (< 300 meV) spin fluctuations is roughly 
1 — 2nh, so one may expect to see relatively little 
change in overall spectral weight compared to the 
undoped antiferromagnet. 

Measurements on Zn-doped Lai.86Sro.i4Cu04 
have shown that the positions and Q-widths of 
the IC peaks are essentially unaffected by the 
presence of the Zn p7[ | . Similarly, doping Zn into 
YBa2Cu306.6 does not change the Q dependence 
of the scattering, but it does shift spectral weight 
down to low-energy, so that x"i^) looks more like 
that of an antiferromagnet [^. This suggests 
that doping Zn into the Cu02 planes destroys su- 
perconductivity by pinning the stripes. 

Recently there has been considerable excite- 
ment over observations of resonant enhancements 
of the magnetic scattering below Tc [19 20 1. If 
the enhanced scattering is assumed to come from 
intinerant electrons, then it is certainly quite re- 
markable. On the other hand, if the magnetic 
scattering all comes from antiferromagnetic do- 
mains between charge stripes, the change in the 
scattering might in some way be associated with a 
reduction in damping of the spin fluctuations due 
to the onset of superconductivity. In any case, 
it will be interesting to see how models of these 
effects evolve. 
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